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Abstract 

We review our knowledge of the most basic properties of the AGN obscuring region - its 
location, scale, symmetry, and mean covering factor - and discuss new evidence on the 
distribution of covering factors in a sample of ~ 9000 quasars with WISE, UKIDSS and 
SDSS photometry. The obscuring regions of AGN may be in some ways more complex 
than we thought - multi-scale, not symmetric, chaotic - and in some ways simpler - with no 
dependence on luminosity, and a covering factor distribution that may be determined by 
the simplest of considerations - e.g. random misalignments. 



1 Introduction 

Since the 1980s, the standard unified scheme for Active Galactic Nuclei (AGN) has included 
an optically and geometrically thick obscuring region, the viewing angle towards which de- 
termines whether we see a Type I or Type II AGN |17[ El [14"1 [37] . Because it must be very 
thick (H/R ~ 1) and is generally assumed to be rotating, the obscuring region has become 
known as the "torus" even though its exact geometry is unclear. Since Krolik and Begelman 
1988 [15] it has been known that such a structure is hard to maintain, and must be clumpy 
or filamentary. Models with a a wedge-like distribution of clumps have had some success 
in explaining the AGN spectral energy distribution (SED) in the MIR (eg [231 E El 134"]). 
However such models are purely phenomenological. What have we learned about the actual 
state of this material? In this conference paper, we look at what the observational evidence 
tells us about the most basic properties of the obscuring region - its location, scale, symmetry, 
and covering factor, and especially the distribution of covering factors. 
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2 Location of the obscuring region 

Much evidence shows that AGN obscuration is actually a multi-scale, multi-component phe- 
nomenon. A good recent summary is in Elvis (2012) [5]. X-ray absorption is often variable 
indicating that it occurs on scales of < 10 3 i?s although there may also be absorption on 
"torus" scales. The characteristic SED peak in the MIR on the other hand suggests optically 
thick dust on scales of ~ W 5 Rs ~ pc x M#/10 8 Af Q - the traditional "torus" component. 
MIR interferometry is just begining to probe these scales ([20, 36J; see also the contribution 
in this volume by Tristram). AGN also seem to be routinely accompanied by star- forming 
discs on a scale of hundreds of parsecs, which can contribute significantly to obscuration (see 
the contributions by Burtscher and by Imanishi) . Finally, the interstellar medium of the host 
galaxy can be crucial, as originally suggested by Keel (1980) [13]. (See also the contribution 
by Goulding) . Although we loosely refer to "the torus" , the true situation may be much more 
complicated than a single neat parsec-scale structure. 

3 Symmetry of the obscuring region 

It is generally assumed that the obscurer is azimuthally symmetric. For a torus supported 
by rotation, or a dust bearing disc wind, this is what we would expect. But if the material 
is part of an incoming streamer (29] or a warped disc [8] it will be not be symmetric. The 
observational sign of this in projection on the sky would be misalignments between various 
AGN structures. Such misalignments in the archetypal AGN NGC 4151, NGC 1068, and Cyg 
A are discussed by |18j . Another example is Cir X-l, a cartoon of which is shown in Fig. [TJ 
Like in NGC1068, the characteristic position angle on the sky seems to change systematically 
with radial scale, a possible sign of incoming material warping to meet the nuclear axis. It 
seems that the "torus" is not necessarily a neat symmetrical doughnut structure. It is more 
likely to be chaotic or systematically distorted and asymmetrical. 

4 Average covering factor 

There are two ways we can estimate the global covering factor of the obscuring region - by 
counting AGN types, or by looking at what fraction of the AGN power is re-processed into 
the IR. The fraction of obscured AGN has been discussed many times, and recently reviewed 
by Lawrence and Elvis (2010) |18j . who find that the fraction of heavily obscured objects 
is 0.58 ± 0.03 and the fraction of lightly obscured objects, where we can still see the BLR, 
is 0.15 ± 0.05. The latter are probably host-obscured objects, so that the fraction of ~ 0.6 
relates to the parsec-scale torus-like region. 

On the other hand, we know that the IR emission constitutes around 30% of the 
bolometric output of the average AGN [30, 6j, suggesting that this is the fraction re-processed 
by the obscuring material. It is interesting that the reprocessed fraction (0.3) agrees roughly 
with the obscured fraction (0.6), but not exactly. We can only measure the re-processed 
fraction for Type I AGN; so this is an indication that the typical properties of Type I and II 
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Figure 1: Cartoon indicating structures in the Circinus galaxy on sub-arcsec scales. A large 
scale radio flow ^ is closely orthogonal to the BLR polarisation 125]/ , indicating the true 
nuclear axis. The inner maser disc fEj matches this, but rotates at large radii. The MIR 
structure [36] can be modelled with two components, with the outer part at a different PA 
from the inner part. On a scale of a few arcsec, there is a one-sided OIII cone. The jet axis 
is inside this cone, but does not bisect it, as shown by JEj. 

AGN are not the same, a point we return to below. 

It is very striking that the obscuring region covers half the sky as seen from the central 
engine. As well as being scarily large, and so hard to explain in rotational models, it is a 
curious number - why not 1%, or 99% ? It is a big clue that we should be looking at either 
disc wind models (eg [3] ) or warped disc models (eg [26J [18] ) both of which naturally give 
rise to covering factors of roughly this order. 



5 Luminosity dependence of covering factor 

There have been many discussions of the possible dependence of obscuration on AGN lu- 
minosity (eg (T71 HU [331 En]; see also the contribution by Toba). There are a number of 
confusing factors and observational problems, such as where to place LINERs and weak-lined 
radio galaxies, the difficulty of getting good classifications in faint and high redshift samples, 
and the conflation of optical obscuration and X-ray absorption, which don't always go to- 
gether. The evidence was reviewed by [18] . who came to a puzzling conclusion. For optical, 
radio, and IR selected samples, there seems to be no effect, with the fraction of Type II AGN 
staying the same at roughly 0.6 over several decades of luminosity. However, X-ray samples 
unambigously show a decline in obscured fraction with increasing luminosity, from 0.8 to 0.1 
over 4-5 decades in X-ray luminosity. One possibility is that the X-ray effect is caused by 
common but unrecognised partial covering by Compton thick material in objects which are 
mistakenly thought to be absorbed but Compton-thin, which would systematically bias the 
apparent luminosities of such objects (Mayo and Lawrence in preparation; see Fig^|. 
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Figure 2: Possible explanation of the luminosity dependence of X-ray obscured fraction, from 
Mayo and Lawrence in preparation. The data are from IlOj . The curves represent various 
models involving populations with varying degrees of partial covering. There is no intrinsic 
variation of obscured fraction with luminosity; the effect in these models is caused by partial 
covering biasing the apparent luminosity. 



6 Probability distribution of covering factor 
6.1 Evidence for a range of covering factors 

The simplest unified scheme assumes that all AGN are identical, with apparent differences 
being due to the accident of viewing angle. The simplest variation assumes that all AGN 
are the same on average, but that there is a distribution of some properties, for example the 
covering factor of obscuring material. As first pointed out in [16], this will lead to Type I and 
II AGN having apparently different properties; objects which happen to have larger covering 
factors, when oriented at random, are more likely to be observed as Type II AGN than 
objects which have lower covering factors. In general if the probability density distribution 
of covering factors C is P(C) then Type II objects will have P{H) = CP(C) and Type I 
objects will have P(I) = (1 - C)P(C) . 

An example is shown in Fig. [3j which shows that for a given radio power, Type I 
AGN have stronger OIII emission than Type II AGN. On the assumption that radio power 
is an indicator of the true UV power, the effect could be produced by a single population 
because objects with a large covering factor will have a smaller open-cone angle, and so 
relatively weaker OIII. The distribution of covering factors can also potentially explain why 
the covering factor implied by the IR luminosity of Type I AGN is half as big as that suggested 
by the Type II/Type I ratio - objects with smaller covering factors are more likely to observed 
as Type I AGN. It could also be relevant to the claim that Type II AGN on average have 
more prominent dust lanes in their circumnuclear regions \21\ [12] . 
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Figure 3: Comparison of radio vs OIII power for Type I and Type II AGN. Taken from 116^ 

6.2 Predictions for the range of covering factors 

A range of covering factors can arise fairly naturally in disc wind models, warped disc models, 
and chaotic accretion models. The only model so far however which has made a concrete 
prediction for the distribution is the simple misaligned disc model |18j . This is only a first 
step towards a proper physical model, but has the advantage of clear predictions with no free 
parameters. The hypothesis is that the direction of incoming material is completely random 
with respect to the nuclear axis, and that the flow warps at some intermediate radius. The 
first simple variant assumes that during the warp the line of nodes rotates fully. This produces 
far too many obscured objects. The second simple variant assumes that there is no rotation of 
the line of nodes; the disc just tilts, with no twist. This predicts fJ-(C a u) = 0.5, /x(C/) = 0.35 
and n{Cn) = 0.7. This very crude model roughly gets correct the overall obscured fraction, 
the typical IR power, and the factor 2 difference in OIII strength. 

6.3 Modeling WISE-UKIDSS-SDSS quasars 

The WISE survey [38] has given us the oppportunity to construct IR-UV SEDS for large 
samples of AGN. In [28J we have used this to study the variation in relative IR to UV 
strength, and hence covering factor, in luminous quasars. We start with the SDSS-DR7 
quasars ([31j|32]; 105,783 objects). We then restrict to > 10 46 erg s _1 and z < 1.5 
(26,927 objects). We then crossmatch with the UKIDSS LAS YJHK survey [IS], which 
covers 4000 sq.deg., giving 9,230 objects - all the SDSS quasars are easily seen in UKIDSS. 
Of these, 9,112 are detected in at least one WISE band, but just 3,831 in all four WISE 
bands. The latter is our prime sample, but we have carefully modelled selection effects by 
use of the larger sample. 

We then fit our SEDs with three empirical components, (i) The mean optical-UV SED 
from [6], representing the accretion disc; (ii) hot dust modelled as a single black body; and 
(iii) a torus component chosen from a library of clumpy torus models [21] . We do not intend 
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Figure 4: Distribution of luminous quasar covering factors, derived from WISE- UKIDSS- 
SDSS SEDs, taken from Roseboom et al [28]. In the upper panel, the black histogram is 
the raw data, and the dashed curve is the version corrected for incompleteness. The blue 
histogram is from the sample of \2T$ . In the lower panel, the red curves are the prediction 
from the tilted disc model fT8\/ ; the dot-dashed line includes the selection correction. 

to extract meaningful physical parameters for the models; they are just an objective way to 
join the dots and characterise the IR luminosity. We then estimate the covering factor as the 
ratio of IR luminosity (both components) to the bolometric luminosity. 

The result is shown in Fig. |4j The observed distribution is approximately a Gaussian 
in logC with [i = —0.41 and a = 0.20, i.e. a mean value fj,c = 0.39 and a one sigma range 
in C of a factor 2.5. This is the distribution in C for Type I AGN; but we can now deduce 
the distribution for all AGN from which this drawn, which has /ic(all) = 0.64, in fairly good 
agreement with the observed obscured fraction, which suggests fic = 0.58. Note that this 
result is for high luminosity quasars, not low-luminosity AGN; it is in strong disagreement 
with the Spitzer-based claim of Treister et al (2008) that the covering factor is ~ 0.04 at 
similar bolometric luminosities. Fig. [4] also shows the prediction from the simple warped disc 
model referred to above. It is in good but not perfect agreement. Note that this prediction 
has no free parameters. It would be good to see predictions from other possible models, such 
as disc wind models. 

6.4 Hot dust behaviour 

In [28], we fitted a hot dust component as well as a cooler torus component. It is clear that 
as well as the overall covering factor having a range of values, the relative strength of the 
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Figure 5: Correlation of the ratio of near-IR and total IR luminosities with various quanti- 
ties: bolometric luminosity, IR luminosity, and overall covering factor. Taken from 128^ . 

hot component also varies considerably. Fig. [5] shows the behaviour of ratio of the l-5//m 
luminosity to the integrated IR luminosity. Two interesting things stand out. The first is that 
the contribution of the hot dust is quite large - typically 30% of the total. This is a significant 
challenge for clumpy torus models and may require a separate component (see \22\ [M] . The 
second interesting thing is that the hot dust fraction anti-correlates with overall covering 
factor - objects with low covering factor have stronger hot dust. We speculate that this may 
be a geometrical effect. The dust nearest the centre will be the hottest and also the most 
likely to be optically thick, so that there may bne a kind of hot inner wall. For geometrically 
thin objects this wall may be preferentially seen face on, whereas for geometrically thick 
objects it may be seen preferentially edge on, as we are looking down the funnel. 



7 Conclusions 



Our main conclusions are as follows, (i) Obscuration is multi-scale and complex, (ii) The 
axisymmetry of the obscuring region is unclear, (hi) The mean covering factor is ~0.6. (iv) 
Evidence on luminosity dependence is conflicting, but it seems most likely that there is no 
intrinsic dependence of covering factor on luminosity, (v) The distribution of covering factors 
is roughly log-Gaussian, with two-thirds within a factor of 2.5. (vi) A simple tilted disc model 
is a good fit to the distribution of covering factors, (viii) Objects with larger covering factor 
have weaker hot dust. 
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